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[']  A  set  of  sensitivity  experiments  are  performed  with  a 
two-way  coupled  and  nested  ocean-atmosphere  forecasting 
system  in  order  to  deconvolve  how  dense  phytoplankton 
stocks  in  a  coastal  embayment  may  impact  thermal  energy 
exchange  processes.  Monterey  Bay  simulations  parameter¬ 
izing  solar  shortwave  transparency  in  the  surface  ocean  as  an 
invariant  oligotrophic  oceanic  water  type  estimate  consis¬ 
tently  colder  sea  surface  temperature  (SST)  than  simulations 
utilizing  more  realistic,  spatially  varying  shortwave  attenu¬ 
ation  terms  based  on  satellite  estimates  of  surface  algal 
pigment  concentration.  These  SST  differences  lead  to  an 
~88%  increase  in  the  cumulative  turbulent  thermal  energy 
transfer  from  the  ocean  to  the  atmosphere  over  the  three 
month  simulation  period.  The  result  is  a  warmer  simulated 
atmospheric  boundary  layer  with  respective  local  air  tem¬ 
perature  differences  approaching  ~2°C.  This  study  suggests 
that  the  retention  of  shortwave  solar  flux  by  ocean  flora  may 
directly  impact  even  short-term  forecasts  of  coastal  meteo¬ 
rological  variables.  Citation:  Jolliff,  J.  K.,  T.  A.  Smith,  C.  N. 
Barron,  S.  deRada,  S.  C.  Anderson,  R.  W.  Gould,  and  R.  A.  Amone 
(2012),  The  impact  of  coastal  phytoplankton  blooms  on  ocean- 
atmospherc  thermal  energy  exchange:  Evidence  from  a  two-way 
coupled  numerical  modeling  system,  Geophys.  Res.  Lett.,  39, 
L24607,  doi:  10.1 029/20 12GL053634. 

1.  Introduction 

[2]  Solar  irradiancc  is  the  principal  source  of  heat  for  the 
surface  ocean.  After  transmission  through  the  intervening 
atmosphere  and  reflection  at  the  sea  surface,  the  penetration 
of  the  remaining  solar  irradiance  is  dependent  upon  the 
absorption  and  scattering  properties  of  both  seawater  and 
the  dissolved  and  suspended  particulate  materials  therein. 
The  bulk  shortwave  solar  irradiant  intensity  encompasses  a 
broad  spectral  range  (~280  to  2500  nm).  Absorption  of  this 
energy  at  infrared  and  longer  wavelengths  (>700  nm)  tends  to 
be  dominated  by  pure  water  [Morel  and  Antoine,  1994;  Smith 
and  Baker,  1981].  However  towards  the  blue  end  of  the 
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ocean  spectral  transparency  window,  phytoplankton  pig¬ 
ments  harvest  the  penetrative  photon  flux  for  photosynthesis. 
Biophysical  rate  limitations  render  this  harvest  remarkably 
inefficient;  most  of  the  energy  is  dissipated  as  heat  Hence  the 
direct  consequence  of  a  variable  surface  ocean  phytoplankton 
biomass  is  the  corresponding  modulation  of  Z^,  the  effective 
depth  over  which  95%  of  the  in- water  penetrative  solar  flux 
is  absorbed  [Morel  and  Antoine,  1994].  Indeed,  global  empir¬ 
ical  relationships  between  surface  chlorophyll-a  concentration 
and  Zfj j  may  be  derived  [see  Morel  et  ai,  2007].  Subsumed  in 
this  simplified  scenario  are  the  associated  and  presumptively 
co-varying  phytoplankton  degradation  materials  that  also 
absorb  light  significantly  in  the  ultraviolet  and  shorter  visible 
wavelengths  [e.g.,  Carder  et  al,  1989]. 

[3]  The  confounding  factors  of  surface  ocean  heat  losses, 
three-dimensional  advection  and  diffusion,  and  variable  solar 
forcing  have  prompted  investigators  to  use  numerical  ocean 
circulation  models  to  investigate  the  manifest  ramifications 
of  the  “ocean  biological  feedback”  (OBF)  effect  outlined 
above.  These  experiments  have  demonstrated  non-negligible 
simulated  SST  sensitivity  to  the  surface  ocean’s  bio-optical 
variability  [e.g.,  Anderson  et  al.,  2007;  Cahill  et  ai,  2008; 
Oschlies,  2004;  Wu  et  al,  2007].  SST,  in  turn,  is  known  to  be 
a  key  variable  that  determines  the  magnitude  and  variability 
of  turbulent  air-sea  transfers  of  heat;  these  exchange  pro¬ 
cesses  then  further  modulate  SST.  What  remains  largely 
unresolved  is  an  understanding  of  how  OBF  may  be  modu¬ 
lating  the  fully  coupled  air-sea  interaction  and  the  time- 
dependent  evolution  of  thermal  variables  in  the  upper  ocean 
as  well  as  the  lower  atmosphere. 

2.  Methods 

[4]  To  begin  to  address  the  fundamental  impact  of  OBF  on 
thermal  exchange  processes  at  the  interface  of  planetary 
boundary  layers,  a  flilly  coupled  ocean-atmosphere  modeling 
system  is  required.  To  that  end,  we  perform  numerical 
experiments  using  the  Coupled  Ocean-Atmosphere  Mesoscale 
Prediction  System  (CO  AMPS®,  a  registered  trademark  of  the 
Naval  Research  Laboratory),  a  nested  numerical  modeling 
system  developed  at  the  Naval  Research  Laboratory  (NRL) 
that  allows  for  a  two-way  exchange  of  information  between 
the  atmospheric  and  oceanic  numerical  forecasting  compo¬ 
nents.  The  COAMPS  configuration  for  our  Monterey  Bay, 
California  study  region  consists  of  two  independent  three- 
dimensional  variational  (3DVAR)  atmosphere  and  ocean 
data  assimilation  systems  and  a  coupler  that  functions  as  a 
router  to  distribute  the  model  forecast  fields  between  the 
oceanic  and  atmospheric  forecast  components.  The  atmo¬ 
spheric  and  oceanic  model  coupling  occurs  via  the  upper- 
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Figure  1.  (a)  COAMPS  nested  model  configuration:  outer  atmospheric  model  nests  begin  at  27  km  horizontal  resolution 
and  nest  iteratively  down  to  9,  3,  and  1  km.  The  inner  ocean  model  nest,  at  500  meter  horizontal  resolution  and  centered 
upon  Monterey  Bay,  California,  is  indicated  by  the  dashed  line,  (b)  Composite  MODIS  image  for  the  surface  chlorophyll-^ 
product;  the  inner  nest  ocean  model  domain  is  indicated  by  the  red  box.  (c)  Bathymetry  map  for  Monterey  Bay,  California 
with  reference  locations  as  listed  in  Table  I  for  Monterey  Bay  Aquarium  Research  Institute  (MBARI)  permanent  moorings 
M 1  and  MO,  and  NRL  temporary  moorings  S 1 02  and  A2  (red  asterisks).  Land  locations  for  California  Irrigation  Management 
Information  System  (CIM1S)  meteorological  stations:  (209)  Watsonville  West  and  (19)  Castroville  are  indicated  with  a 
square.  The  Elkhom  Slough  (Elk)  reference  site  is  approximate  to  the  National  Estuarine  Research  Reserve  System  (NERRS) 
weather  station. 


most  oceanic  model  grid  cell  temperature  and  the  lowest  grid 
cell  atmospheric  model  variables  (temperature,  humidity, 
wind  velocity,  pressure,  and  radiative  fluxes).  Bulk  fluxes 
of  heat  energy  exchange  are  calculated  following  the  Cou¬ 
pled  Ocean-Atmosphere  Response,  version  3  (COARE  3.0) 
scheme  [Fairall  et  al,  1996]. 

[5]  The  Navy  Coastal  Ocean  Model  (NCOM)  [Barron  et  al. , 
2006]  serves  as  the  ocean  model  component  with  assimilation 
performed  via  the  Navy  Coupled  Ocean  Data  Assimilation 
(NCODA)  3DVAR  system  [Cummings,  2005].  The  COAMPS 
domain  configuration  for  Monterey  Bay  (Figure  la)  consists 
of  a  quadruple-nested  atmospheric  domain  with  27,  9, 3,  and 
1  km  horizontal  spacing.  The  NCOM  configuration  consists 
of  an  inner  nest  with  500  meter  horizontal  grid  spacing.  It  is 
initialized  using  fields  interpolated  from  the  global  (1/8°) 
NCOM  results  [Kara  et  al ,  2006];  similar  fields  serve  as 
boundary  conditions  during  the  4-month  simulation.  This 
inner  nest  is  executed  in  data-assimilative  mode  during  a 
1  month  spin -up  period  (beginning  on  I  April  2008)  and 
then  in  non-data-assimilative  mode  (or  ‘free-run*  mode) 
for  the  remaining  three  months  of  the  simulation  (ending 
31  July  2008). 

[6]  Verification  and  validation  of  the  COAMPS  forecast¬ 
ing  system  may  be  found  elsewhere  [Doyle  et  al. ,  2009;  Small 
et  al ,  2012];  here  we  focus  on  the  modeling  system’s  sensi¬ 
tivity  to  changes  in  the  oceanic  shortwave  optical  attenuation 
that  arise  as  a  result  of  surface  ocean  phytoplankton  vari¬ 
ability.  Previous  COAMPS  applications  in  the  Monterey  Bay 
region  have  suggested  that  a  40%  reduction  factor  for  surface 
shortwave  irradiance  is  warranted  [Shulman  et  al ,  2007]. 
Such  a  reduction  compensates  for  the  inability  of  COAMPS 
to  simulate  the  specific  reflective  effects  of  low-lying  marine 
stratus  clouds  [Doyle  et  al ,  2009].  Here  we  present  results 
from  a  ‘reduced  shortwave’  series  of  numerical  experiments. 


[7]  Specific  to  the  simulated  attenuation  of  solar  short¬ 
wave  in  the  ocean  model,  the  default  parameterization  for 
COAMPS  is  the  “Jerlov”  [Jerlov,  1976]  oceanic  water  Type 
1A  implemented  in  a  two-band  shortwave  extinction  scheme 
popularized  for  physical  ocean  models  by  Paulson  and 
Simpson  [1977]  (herein,  ‘PS77-LA’).  The  main  feature  to 
be  distinguished  here  is  that  PS77-LA  employs  time/spaee 
invariant  scalar  quantities  associated  with  an  ideal  oligotro- 
phic  state,  thereby  ignoring  any  potential  optical  variability 
in  surface  waters  over  the  model  domain.  The  PS77-1A 
parameterization  is  used  for  the  control  model  simulation 
case  (SI). 

[s]  An  upgraded  representation  of  shortwave  attenuation 
is  examined  in  the  alternative  simulation  ease  (S2),  which 
replaces  the  ideal  oligotrophic  conditions  of  PS77-LA  with 
a  satellite-based  attenuation  scheme  that  allows  recent  obser¬ 
vations  to  guide  attenuation  levels.  For  this  purpose,  a  4-month 
composite  is  prepared  based  on  sea  surface  chlorophyll 
(mg  m~3)  estimates  derived  from  the  Moderate  Resolution 
Imaging  Spcctrorad iometer  (MODIS)  sensor  and  mapped  to 
the  inner  ocean  model  nest  (Figure  lb);  each  pixel  in  the 
composite  product  averages  all  clear  observations  of  that  pixel 
over  the  4-month  period.  The  two-dimensional  surface  chlo¬ 
rophyll  fields  are  read  into  NCOM  for  calculations  of  the 
depth-dependent  transmission  of  shortwave  irradiance  using 
the  Lee  et  al  [2005]  method.  This  scheme  divides  total  solar 
shortwave  into  Photosynthetically  Available  Radiation  (PAR; 
350-700  nm)  and  longer  (>700  nm)  spectral  components.  The 
depth-dependent  attenuation  coefficient  for  the  integrated 
PAR  spectral  range  [KPAR(z)]  is  then  calculated  based  on  the 
surface  absorption  and  scattering  coefficients.  These  optical 
coefficients  are  inferred  from  the  satellite  surface  chlorophyll 
product  using  algorithms  based  on  published  ‘Case  F  rela¬ 
tionships  [Morel,  1988;  Morel  et  al,  2007].  A  separate 
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Figure  2.  (a)  Simulated  shortwave  radiation  transmittance  using  the  PS77-IA  parameterization  (SI)  and  the  satellite-based 
method  described  in  the  text  (S2)  for  the  MO  reference  location  in  Monterey  Bay.  (b)  The  temperature  profile  bias  for  each 
month  at  the  MO  location  is  calculated  by  subtracting  the  mean  monthly  temperature  profile  of  SI  from  the  mean  profile  of 
S2,  the  May  and  June  profiles  are  nearly  identical,  the  July  profile  is  labeled,  (c)  The  3-hourly  SST  results  for  S2  (red)  and  SI 
(black),  (d)  The  3-hourly  SST  difference  (S2  -  SI),  the  smoothed  trend  line  (solid  black)  and  the  bias  (dashed  red)  are  shown, 
(e)  The  ratio  of  SI  depth- integrated  thermal  energy  content  (for  the  upper  10.45  meters)  to  the  same  as  calculated  for  S2.  As 
in  Figure  2d,  the  smoothed  series  are  shown  as  well  as  the  respective  biases  (differences  of  the  respective  time  series  means; 
dashed  red  line). 


analysis  of  in  situ  optical  data  [Jolliff  et  aly  2012]  confirmed 
these  relationships  to  be  reasonable  for  the  Monterey  Bay 
region. 

[9]  The  nested  modeling  system  domain  is  centered  upon 
the  Monterey  Bay  region  of  the  U.S.  West  Coast  (Figure  la). 
Monterey  Bay  is  a  --40  km  wide  open  embayment  known  for 
very  high  rates  of  biological  productivity  [Wilkerson  et  al.y 
2000],  particularly  during  the  spring  to  summer  upwelling 
season  [Chavez  et  a/.,  1991].  As  is  typical  of  eastern  bound¬ 
ary  upwelling  systems,  colder  subsurface  waters  emerge  at 
the  surface  along  the  coastal  divergence,  and  rapid  warming 
of  these  waters  is  concomitant  with  biological  utilization  of 
dissolved  nutrients  leading  to  persistent  coastal  eutrophica¬ 
tion  (defined  here  as  >2.0  mg  m“3  surface  chlorophyll-a) 
[Kudela  and Dugdale,  2000].  Simulated  land/sea  temperature 
comparison  reference  points  arc  selected  from  the  coupled 


model  domain  (Figure  lc)  in  order  to  assess  the  thermal 
impact  of  phytoplankton  blooms  inside  Monterey  Bay.  These 
reference  points  correspond  to  temporary  and  permanent 
marine  mooring  locations  and  land  surface  meteorological 
observation  stations. 

3.  Results 

[10]  In  comparison  to  the  satellite-based  shortwave  atten¬ 
uation  scheme,  the  PS77-IA  parameterization  overestimates 
Zh,  by  33  meters  at  mooring  M0  (Figure  2a).  The  important 
distinction  between  the  SI  and  S2  attenuation  profiles  is  that 
under  the  satellite-determined  chlorophyll  concentrations  of 
S2,  the  upper  decameter  of  the  water  column  absorbs  most  of 
the  penetrative  shortwave.  The  result  of  this  curtailed  pene¬ 
trative  flux  is  evident  in  the  temperature  profile  bias,  i.e., 
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Table  1.  SI  and  S2  Comparison  Statistics* 


SI  versus  S2  for  May-Junc-July 

Ml 

M0 

NRL-S102 

NRL- A2 

Latent  Heat  Flux  Bias 

3.79  Wm  2 

4.7  Wm  2 

3.9  Wm"2 

3.46  W  m"2 

Lalcnt  Heat  Flux  Amplitude,  r2 

1.13,0.88 

1.24,0.88 

1.27,0.88 

1.14,0.89 

Total  Latent  Increase  (%) 

30.8 

40.4 

34.3 

32.6 

Turbulent  Heat  Flux  Bias 

4.33  Wm  2 

5.22  Wm  2 

4.11  Wm  2 

3.9Wm'2 

Turbulent  Heat  Flux  Amplitude,  r2 

1.12,0.82 

1.27,0.83 

1.29,0.88 

1.11,0.85 

Total  Turbulent  Heat  Flux  Increase  (%) 

66.7 

63.1 

44.2 

87.9 

Air  Temp.  (2  m)  Bias 

0.39°C 

0.5I°C 

0.60°C 

0.39°C 

Air  Temp.  (2  m)  Amplitude,  r2 

1.07,0.98 

1.08,0.97 

1.12,0.98 

1.10,0.98 

Maximum  Air  Temp.  (2  m)  Diff. 

1.5°C 

1.6°C 

1.6°C 

2.2°C 

SST  (1  m)  Bias 

0.45°C 

0.57°C 

0.64°C 

0.43°C 

SST  ( 1  m)  Amplitude,  r2 

1.05,0.96 

1.09,0.96 

1.13,0.97 

1.10,0.96 

Maximum  SST  (1  m)  Difference 

1.55°C 

I.70°C 

1.77°C 

1.60°C 

*For  each  of  the  (4)  marine  reference  points  indicated  in  the  lop  row,  time  series  of  simulated  thermal  variables  (hourly  results  for  May,  June,  and  July) 
arc  extracted  from  each  respective  simulation  (SI  and  S2).  Bias  refers  to  the  difference  between  the  means  calculated  for  each  simulation  time  series:  mean 
(S2  results)  -  mcan(Sl  results).  Amplitude  is  the  ratio  of  standard  deviations:  SD(S2ySD(Sl).  The  F  statistic  is  the  square  of  the  linear  correlation 
coefficient  between  the  two  time  series.  The  percentage  increases  (%)  in  the  heal  flux  terms  arc  calculated  from  the  total  lhcrmal  energy  passing  from 
the  ocean  to  the  atmosphere  (F)  during  the  respective  three  month  free-run  simulation  periods  for  SI  and  S2  following  the  formula:  [F(S2)  —  F(S1)J/ 
F(S1)  *  100.  The  turbulent  heat  flux  is  the  sum  of  the  latent  and  sensible  heat  fluxes. 


the  difference  between  the  monthly  average  temperature 
profiles  (Figure  2b).  The  July  monthly  bias  is  positive  (S2  is 
warmer)  by  0.7°C  in  the  upper  decamater,  below  this  depth 
S2  is  colder  by  ~0.I-0.2°C.  This  comparative  redistribution 
of  thermal  energy  is  to  be  expected  based  on  the  different 
shortwave  attenuation  profiles:  neglecting  advcction  and 
vertical  mixing,  the  direct  solar  heating  rate  is  inversely 
proportional  to  the  heating  layer  thickness,  Zh|. 

[u]  The  simulated  temporal  SST  trends  at  reference 
mooring  MO  are  indicative  of  the  spring  transition  from  a 
well-mixed  water  column  to  increasing  thermal  stratification 
over  the  three  months  of  the  free-run  simulation  period 
(May-July;  Figure  2c).  SSTs  in  both  simulations  increase  by 
at  least  ~4°C  over  the  transition  period.  Nonetheless,  the 
three-month  S2  warm  bias  is  comparatively  consistent  at 
+0.57°C  with  the  instantaneous  SST  difference  ranging  as 
high  as  ~  1 .7°C  (Figure  2d).  The  upper  decameter  integrated 
thermal  energy  in  SI  fell  as  much  as  ~10%  lower  than  in  S2 
(Figure  2c)  with  a  mean  difference  of  3%.  The  comparatively 
consistent  SST  bias  over  the  free-run  simulation  period  in 
conjunction  with  the  high  correlation  of  the  SI  and  S2  SST 
results  (r2  =  0.96)  indicates  that  the  main  impact  of  the  OBF 
effect  in  the  S2  simulation  is  a  change  in  the  thermal  exchange 
of  energy  between  the  simulated  ocean  and  atmosphere. 

[12]  Indeed,  the  cumulative  total  turbulent  thermal  energy 
transfer  (quantified  as  the  sum  of  the  latent  and  sensible  heat 
flux  to  the  atmosphere)  increases  by  87.9%  at  the  model  loca¬ 
tion  corresponding  to  NRL  temporary  mooring  A2  (Figure  lc 
and  Table  1).  Similarly,  the  corresponding  increase  for  MO  is 
63%.  Wind  speeds  arc  largely  the  same  in  the  two  simulations, 
although  some  fine- scale  local  mismatches  are  evident  (wind 
results  not  shown).  This  indicates  that  the  increased  thermal 
energy  in  the  surface  layer  of  S2  drives  simulated  differences 
in  latent  transfer  via  an  increased  SST  and  changes  in  the 
turbulent  transfer  coefficients  of  the  Fair  all  et  al  [  1 996]  bulk 
flux  scheme.  Accordingly,  the  lower  marine  air  temperatures 
at  the  mooring  locations  are  elevated  in  the  S2  simulation 
(Table  1).  Both  the  SST  and  lower  air  temperature  biases  and 
maximum  instantaneous  differences  closely  match  one 
another.  The  high  correlation  values  (r2)  indicate  that  poten¬ 
tial  differences  in  advection  (surface  winds,  currents)  are  not 
significantly  impacting  the  simulation  comparisons  at  fixed 
reference  points  inside  Monterey  Bay.  The  slight  increase  in 


amplitude  (the  ratio  of  time  scries  standard  deviations)  is 
indicative  of  a  stronger  diurnal  warming  signal  in  S2. 

[13]  The  manifest  result  of  the  increased  thermal  energy 
transfer  in  S2  is  an  increase  in  the  simulated  air  temperature 
for  the  marine  atmospheric  boundary  layer  (MABL)  over 
Monterey  Bay  and  adjacent  coastal  areas.  Monthly  air  tem¬ 
perature  profile  biases  prepared  for  the  coastal  land  reference 
sites  (Figure  1  c)  indicate  a  peak  bias  for  July  of  +0.5°C  at  the 
Castrovillc  meteorological  reference  station  (Figure  3).  All 
reference  stations  indicate  the  simulated  temperature  biases 
propagate  to  >400  m  height.  More  detailed  examination  of 
the  lower  air  temperature  comparisons  (Figure  3,  bottom) 
reveal  an  enhanced  diurnal  warming  signal  that  exceeds  1°C 
difference  at  the  Castroville  station  on  July  18.  July  daily 
maximum  temperature  differences  are  generally  in  the  range 
of  0.5-1 .0°C  for  the  other  land  stations  (data  not  shown). 

[u]  Hence  over  the  course  of  the  three  month  free-run 
simulations,  the  S2  version  reveals  an  increasingly  warm 
MABL  over  Monterey  Bay  (Figure  4).  This  warm  pertuibation 
is  advected  by  the  prevailing  northwesterly  winds  towards  the 
Salinas  Valley  (Figure  4c).  The  air  temperature  warming  bias 
is  concentrated  inside  Monterey  Bay  and  is  consistent  with  the 
MOD1S  chlorophyll  composite  imagery  indicating  phyto¬ 
plankton  blooms  along  the  interior  coastline  of  Monterey  Bay 
(Figure  lb).  In  the  July  bias  field  (Figure  4c)  an  additional 
warm  bias  center  is  apparent  outside  of  Monterey  Bay  along 
36.8°N  latitude.  This  appears  to  result  from  differences 
between  the  two  simulations  in  the  intensity  of  upwelling 
along  the  coast  north  of  Monterey  Bay.  Further  analysis  of  the 
ocean  circulation  changes  are  beyond  the  scope  of  this  paper, 
but  we  note  that  they  become  increasingly  evident  as  the  free- 
run  ocean  model  solution  moves  farther  away  in  time  from  the 
constraints  of  the  ocean  data- assimilation  period  (April  2008). 

4.  Summary  and  Conclusions 

[15]  Here  we  show  that  simulated  transfers  of  thermal 
energy  between  the  ocean  and  atmosphere  are  significantly 
impacted  by  the  enhanced  optical  attenuation  attendant  to 
coastal  phytoplankton  blooms.  Mechanistically,  the  addi¬ 
tional  solar  shortwave  retained  by  surface  phytoplankton 
stocks  remains  available  to  flux  back  into  the  lower  atmo¬ 
sphere  via  turbulent  heat  flux  losses,  thereby  resulting  in  a 
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Figure  3.  (top)  For  land  locations,  the  mean  air  temperature  profile  from  the  atmospheric  model  simulation  (S2)  is  sub¬ 
tracted  from  the  same  for  SI.  Positive  values  indicate  a  warming  bias  for  S2.  The  reference  station  locations  correspond 
to  those  mapped  in  Figure  le.  (bottom)  The  hourly  lower  air  temperature  difference  (S2  —  SI)  is  shown  for  reference  land 
location  19  (CIMIS,  Castroville).  The  3-month  time  series  bias  is  the  solid  black  line,  the  monthly  biases  (May,  June,  and 
July)  are  indicated  by  the  dashed  red  line. 


comparative  surplus  of  thermal  energy  in  the  MABL — 
energy  that  would  otherwise  be  sequestered  in  the  ocean. 
This  ‘reflux’  of  thermal  energy  then  appears  to  have  a  sig¬ 
nificant  impact  on  simulated  lower  air  temperatures:  up  to 
~1°C  locally  over  land  and  ~2°C  over  water.  These  thermal 
atmospheric  perturbations  are  similar  to  those  obtained  by 
Shell  el  al.  [2003]  using  a  coarse  resolution  coupled  atmo- 
sphere-oeean  general  circulation  model;  however,  here  we 
show  that  this  impact  is  pertinent  to  local  scale  air-sea  fluxes 
on  much  shorter  timescales. 

[i6]  Accordingly,  increasingly  capable  representations  of 
optical  attenuation  due  to  dissolved  and  particulate  seawater 
constituents  appear  justified  for  improving  the  potential 
accuracy  of  coupled  ocean-atmosphere  forecasting  systems. 
In  our  numerical  experiments,  the  surface  chlorophyll-^ 


satellite  product  is  used  to  estimate  the  in  situ  absorption 
coefficients  required  for  our  shortwave  attenuation  compu¬ 
tations.  This  approach  is  reasonable  for  the  Monterey  Bay 
region  since  the  measured  surface  water  total  absorption 
coefficients  in  the  visible  spectral  range  and  the  surface 
pigment  concentrations  are  highly  correlated  in  this  system 
[Jolliff  el  al .,  2012].  For  many  other  coastal  environments, 
however,  the  additional  optical  property  contributions  from 
non-living  organic  matter  as  well  as  suspended  sediments 
may  also  need  to  be  considered  explicitly  and  independently 
of  the  surface  phytoplankton  pigment  distribution.  More 
broadly  applicable  coastal  numerical  modeling  efforts 
should  account  for  these  additional  optical  complexities. 

[17]  Our  model  results  nonetheless  leave  open  the  ques¬ 
tion  of  how  anthropogenic  coastal  eutrophication  may  also 
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Figure  4.  The  spatial  distribution  of  the  respective  lower  air 
temperature  (2  meter)  biases  [mean  (S2)  -  mean(Sl)]  are 
shown  for  (A)  May,  (B)  June,  and  (C)  July  of  the  simulations. 


impact  temperatures  at  the  air-sea  interface.  The  degree 
of  thermal  impact  must  be  gauged  with  respect  to  a  control 
point,  which  in  our  numerical  experiments  is  the  PS77-IA 
two-band  abbreviation  of  the  Jerlov  [1976]  Type  IA  oceanic 
surface  waler.  Our  calculations  suggest  that  the  PS77-LA 
attenuation  scheme  corresponds  to  approximately  0.13  mg 
m~3  surface  chlorophyll,  i.e.,  oligotrophic  marine  conditions. 
Given  such  clear  natural  waters  as  a  control  point,  a  signifi¬ 
cant  increase  in  the  surface  absorption  coefficients  due  to 
phytoplankton  blooms,  significant  accumulation  or  discharge 
of  colored  dissolved  organic  matter  (CDOM),  or  both  events 
simultaneously  may  influence  thermal  variables  at  the  air-sea 
interface.  The  effect  may  be  significant  enough  to  discemibly 
impact  daily  land  surface  temperature  readings  near  the  coast 
in  the  M).5-1.0°C  range. 
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